Spin Dynamics and Orbital State in LaTiOs 
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A neutron scattering study of the Mott-Hubbard insulator LaTiOs (Tn = 132 K) reveals a spin wave 
spectrum that is well described by a nearest-neighbor superexchange constant J = 15.5 meV and a 
small Dzyaloshinskii-Moriya interaction (D = 1.1 meV). The nearly isotropic spin wave spectrum 
is surprising in view of the absence of a static Jahn- Teller distortion that could quench the orbital 
angular momentum, and it may indicate strong orbital fluctuations. A resonant x-ray scattering 
study has uncovered no evidence of orbital order in LaTiOs. 



In the layered cuprates exemplified by the series 
La2- a; Sr 2 ; CUO4+5, the transition from a 3c? 9 antiferro- 
magnetic (AF) insulator at x = S = into an uncon- 
ventional metallic and superconducting state with in- 
creasing hole concentration (x, 6 > 0) has received an 
enormous amount of attention. The magnetic spectra of 
these materials, revealed by inelastic neutron scattering, 
have played a key role in efforts to arrive at a theoreti- 
cal explanation of this transition. The pseudocubic per- 
ovskite Lai-zSr^TiOs+i undergoes an analogous transi- 
tion from a 3d 1 AF insulator at x = S = to a metallic 
state with increasing hole concentration Q. In the ti- 
tanates, however, the metallic state shows conventional 
Fermi liquid behavior, and no superconductivity is found 
jij. Momentum- resolved probes such as angle-resolved 
photoemission spectroscopy and inelastic neutron scat- 
tering have thus far not been applied to the titanates, and 
the origin of the very different behavior of the metallic 
cuprates and titanates is still largely unexplored. Here 
we report an inelastic neutron scattering and anomalous 
x-ray scattering study of the parent compound of the 
titanate series, LaTiC>3, that provides insight into the 
microscopic interactions underlying this behavior. 

Orbital degrees of freedom, quenched in the layered 
cuprates by a large Jahn- Teller (JT) distortion of the 
CuC>6 octahedra, are likely to be a key factor in the phe- 
nomenology of the titanates. While the TiOe octahedra 
are tilted in a GdFeC>3-type structure, their distortion 
is small and essentially undetectable in neutron powder 
diffraction experiments on LaTiC>3 (Ref. |2j). The crys- 
tal field acting on the Ti 3+ ion is therefore nearly cu- 
bic, and heuristically one expects a quadruply degenerate 
single-ion ground state with unquenched orbital angular 
momentum opposite to the spin angular momentum due 
to the spin-orbit interaction. In other perovskites such 
as LaMnC>3, such spin-orbital degeneracies are broken 
by successive orbital and magnetic ordering transitions 



|3[. In the orbitally and magnetically ordered state of 
LaMnC>3, the spin wave spectrum is highly anisotropic 
reflecting the different relative orientations of the or- 
bitals on nearest-neighbor Mn atoms in different crys- 
tallographic directions Q. 

The reduced ordered moment (/io ~ 0.45^b, Ref. jql) 
in the G-Type AF structure of LaTi0 3 (inset in Fig. y]) 
at first sight appears consistent with a conventional sce- 
nario in which the orbital occupancies at every site are 
established at some high temperature, and the magnetic 
degrees of freedom (coupled spin and orbital angular mo- 
menta) order at a lower temperature. Full theoretical 
calculations, however, generally predict a ferromagnetic 
spin structure for LaTiOs ||. As in LaMn03, the spin 
dynamics of LaTiOs are highly sensitive to the orbital oc- 
cupancies and can provide important information in this 
regard. We find that the exchange anisotropy is small 
and hence inconsistent with the presence of an apprecia- 
ble unquenched orbital moment. At the same time, syn- 
chrotron x-ray scattering experiments have not revealed 
any evidence of reflections showing a resonant enhance- 
ment at the Ti K-edge, unlike other perovskites in which 
orbital order (OO) is present. These observations, along 
with previously puzzling Raman scattering data , indi- 
cate strong fluctuations in the orbital sector of LaTiOs. 

The neutron scattering experiments were conducted on 
the BT2 and BT4 triple axis spectrometers at the NIST 
research reactor and at the IN8 spectrometer at the In- 
stitut Laue-Langevin. For excitation energies up to 25 
meV, we used high resolution configurations with verti- 
cally focusing pyrolytic graphite (PG) (002) monochro- 
mator and PG (002) analyser crystals set for final neu- 
tron energies of 14.7 meV or 30.5 meV and horizontally 
collimated beams at both NIST and ILL. For excitation 
energies of 20 meV and higher, we used a double-focusing 
analyser on IN8 with open collimations, and with a Cu 
(111) monochromator and a PG (002) analyser set for 
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a final neutron energy of 35 meV. PG filters were used 
in the scattered beam to reduce higher order contamina- 
tion. Data obtained on the different spectrometers and 
with different configurations were in good agreement. 

The sample was a single crystal of volume 0.2 cm 3 and 
mosaicity 0.5° grown by the floating zone technique. It is 
semiconducting, and neutron diffraction (Fig. |l|) shows a 
sharp, second-order Neel transition at Tn = 132 K, im- 
plying a highly homogeneous oxygen content (5 ~ 0.01 in 
LaTi0 3+(5 ) Q. The diffraction patter n is consistent with 
the G-type structure found previously and a small 
uncompensated moment (~ 10~ 2 /xb per Ti spin at low 
temperatures) appears below Tn due to spin canting, as 
observed by magnetization measurements. The nuclear 
structure of LaTi03 is orthorhombic (space group Pnma 
H), but the crystal is fully twinned. Because of the 
isotropy of the spin wave dispersions (see below), twin- 
ning did not influence the neutron measurements. For 
simplicity we express the wave vectors in the pseudocubic 
notation with lattice constant a ~ 3.95A. In this nota- 
tion, AF Bragg reflections are located at (h/2,k/2,l/2) 
with h, k, I odd. Data were taken with the crystal in two 
different orientations in which wave vectors of the form 
(h, h, I) or (h, k, (h + k)/2), respectively, were accessible. 

Fig. H shows inelastic neutron scattering data obtained 
in constant-q mode with high resolution near the AF zone 
center. The dispersing peak shown disappears above the 
Neel temperature, thus clearly identifying itself as a spin 
wave excitation. In Fig. ^, constant-energy scans ob- 
tained in the high-intensity configuration with relaxed 
resolution are presented. The profile shapes are strongly 
influenced by the spectrometer resolution, and a decon- 
volution is required to accurately extract the positions 
of the spin wave peaks. For the high-resolution config- 
uration with only vertical focusing we used the stan- 
dard Cooper-Nathans procedure while a Monte-Carlo 
ray-tracing routine was used for the doubly focused ge- 
ometry An analytical approximation to the Ti 3+ 
form factor and the standard intensity factors for AF 
magnons were incorporated in the programs pi. The 
peak positions thus obtained are shown in Fig. 13 in the 
(111) direction. 

A very good global fit to all data was obtained 
by convoluting the resolution function with a sin- 
gle spin wave branch of the generic form ho; = 
zSJ\J (1 + e) 2 — 7 2 where ha; is the spin wave energy, 
z = 6 is the coordination number, S = 1/2 is the 
Ti spin, J = 15.5 ± 1 meV is the isotropic (Heisen- 
berg) part of the nearest-neighbor superexchange p0[ , 
7 = |[cos(g x a) + cos(q y a) + cos(<? 2 a)] with the magnon 
wave vector q measured from the magnetic zone center, 
and the zone center gap is A ~ zSJ\[2t = 3.3 ±0.3 meV. 
The solid lines in Figs. ^-|] result from this global fit and 
obviously provide a good description of all data. Inclu- 
sion of further-neighbor interactions, damping parame- 
ters above the instrumental resolution, or other (nondc- 



generate) spin wave branches did not improve the fit. 

The Heisenberg exchange constant J is in fair agree- 
ment with predictions based on a comparison of the 
Neel temperature with numerical simulations (Tn = 
0.946J//cs ~ 170 K for spins-1/2 on a simple cubic 
lattice p]]). In general, the spin wave gap is deter- 
mined by symmetric and antisymmetric (Dzyaloshinskii- 
Moriya) anisotropy terms in the superexchange matrix, 
by terms originating from direct exchange, by dipolar in- 
teractions, and by the single ion anisotropy. The latter 
two effects are negligible and nonexistent, respectively, in 
spin-1/2 systems. In the GdFe03 structure, antisymmet- 
ric exchange is allowed by symmetry, and its magnitude 
is expected to scale with the tilt angle of the TiOg octahe- 
dra. Because of the large tilt angle (11.5°), we expect this 
effect to dominate over the more subtle direct exchange 
terms. Theories of superexchange anisotropies were 
recently reexamined 13 1 in the light of neutron scattering 
data on the layered cuprates. It was shown that the sym- 
metric and antisymmetric terms are related by a hidden 
symmetry so that the two zone-center spin wave gaps de- 
pend on the relationship between Dzyaloshinkii-Moriya 
(DM) vectors centered on each magnetic bond. Specifi- 
cally, for two-dimensional (2D) spin structures (such as 
the one of La2Cu04) the gaps are degenerate if all DM 
vectors have the same magnitude, and the degenerate 
gaps are nonzero if, in addition, not all vectors have the 
same orientation. The bond-dependent DM vectors for 
LaMnOa (isostructural to LaTiOs) were given in Ref. 
Jl4j]. Although a detailed analysis of the spin dynam- 
ics for these 3D systems has not been reported, we note 
that both of the above criteria are fulfilled for LaTi03. 
In particular, the DM vectors centered on the six Ti- 
O-Ti bonds have different orientations but are expected 
to have the same magnitude because of the practically 
equal bond lengths and bond angles [0. Our observa- 
tion of degenerate but nonzero spin wave gaps therefore 
provides support for the predictions of Ref. |13| in a 3D, 
non-cuprate system. 

Carrying the analogy to the 2D cuprates one step fur- 
ther, we expect that the anisotropy gap is A ~ zSD, 
and hence D ~ 1.1 meV is the net DM interaction per Ti 
spin. Due to spin canting, the net ferromagnetic moment 
per spin should therefore be ~ fioD/2J — 1.5 x 10~ 2 /is, 
in good agreement with the observed value. This sup- 
ports our assumption that the DM interaction provides 
the dominant contribution to the spin wave gap. A mi- 
croscopic calculation of D would be a further interesting 
test of the formalism developed in Ref. [O. 

The small easy-axis anisotropy in the exchange Hamil- 
tonian of LaTi03 is difficult to understand based on sim- 
ple crystal-field considerations for the Ti 3+ ion. The anti- 
symmetric exchange is generally of order (Ag/g)J where 
g is the free-electron Lande factor and Ag is its shift in 
the crystalline environment p"2| . Based on our data, we 
therefore estimate Ag/g ~ 0.05. On the other hand, in 
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the absence of any appreciable static JT distortion as ob- 
served experimentally [|| , one expects that the spin-orbit 
interaction ( A ~ +20 meV splits the t2 g multiplet 
of the cubic crystal field Hamiltonian into a quadruply 
degenerate ground state and a higher-lying Kramers dou- 
blet. In a simple crystal field model, this ground state 
is characterized by an unqucnchcd orbital moment equal 
and antiparallel to the spin moment (Ag/g = 1). More 
elaborate Hartree-Fock calculations ||] do not change this 
picture qualitatively: Even if a static JT distortion at the 
limits of the experimental error bars [0 is included, the 
orbital contribution to the moment remains comparable 
to the spin moment (Ag/g ~ 0.5). There is thus an 
order-of-magnitude discrepancy between the predictions 
of conventional models and the neutron scattering ober- 
vations. The smallness of the spin anisotropy in LaTiC>3 
is underscored by a different comparison: The D/J ra- 
tio of LaTiC>3 differs only by a factor of 3 from that 
of La2CuC>4, whose low-temperature ordered moment is 
in near-perfect agreement with the spin-only prediction. 
Since D/J scales with the tilt angle of the octahcdra 
which is a factor of ~ 3 larger in LaTiC>3, the relative 
magnitude of this quantity in the two materials can be 
accounted for without invoking a large orbital moment 
in LaTiC>3. 

Interestingly, the large discrepancy between the pre- 
dictions of conventional models and the neutron scat- 
tering observations has a close analogy in the electron 
spin resonance (ESR) literature. The description of ESR 
data on Ti 3+ impurities embedded into perovskite lat- 
tices in fact commonly requires g-factors that are much 
more isotropic than predicted by simple crystal field cal- 
culations jl^]. According to a widely used model |l6| , 
this is attributed to the dynamical JT effect where the 
orbital degeneracy is lifted by coupling to zero-point lat- 
tice vibrations. 

While the dynamical JT effect is well established in 
impurity systems, it has thus far not been reported in 
lattice systems which commonly exhibit static, coopera- 
tive JT distortions associated with 00. Anomalous x-ray 
diffraction with photon energies near an absorption edge 
of the transition metal ion has recently been established 
as a direct probe of 00 in perovskites || whose sensitiv- 
ity far exceeds conventional diffraction techniques that 
probe 00 indirectly through associated lattice distor- 
tions. We have therefore carried out an extensive search 
for reflections characteristic of orbital ordering near the 
Ti K-edge (4.966 keV) at beamline X22C at the National 
Synchrotron Light Source, with energy resolution ~ 5 
eV. The experiments were performed at low temperature 
(T=l OK) on a polished (1,1,0) surface of the same crystal 
that was also used for the neutron measurements. No ev- 
idence for resonant reflections at several high symmetry 
positions (such as (^,^,0), the ordering wave vector ex- 
pected for t2 g orbitals with a G-type spin structure |17|) 
was found under the same conditions that enabled their 



positive identification in LaMn03 ||, YTi03 |D| and re- 
lated materials. If 00 is present in LaTi03, we can hence 
conclude that its order parameter is much smaller than in 
comparable perovskites. On general grounds, the reduc- 
tion of the order parameter should be accompanied by en- 
hanced orbital zero-point fluctuations. These may have 
already been detected (though not identified as such): A 
large electronic background and pronounced Fano-type 
phonon anomalies were observed by Raman scattering 
in nominally stoichiometric, insulating titanates and are 
most pronounced in LaTiOs |?J . In the light of our obser- 
vations, it is of course important to revisit these exper- 
iments and rule out any possible role of residual oxygen 
defects, inhomogeneity, etc. The presently available data 
are, however, naturally interpreted as arising from orbital 
fluctuations coupled to lattice vibrations. 

An observation not explained by these qualitative con- 
siderations is the small ordered moment in the AF state. 
If the orbital moment is indeed largely quenched, one 
would naively expect a spin-only moment of 0.85/1 b [fL9| , 
in contrast to the experimental observation of 0.45/is. 
A recently proposed full theory of the interplay between 
the orbital and spin dynamics in LaTiOs has yielded a 
prediction of the ordered moment that is in quantitative 
agreement with experiment ]20| ]. 

In conclusion, several lines of evidence from neutron, 
x-ray and Raman scattering can be self-consistently in- 
terpreted in terms of an unusual many body state with 
AF long range order but strong orbital fluctuations. This 
should be an interesting subject of theoretical research. 
The orbital fluctuations are expected to be enhanced in 
the presence of itinerant charge carriers and therefore 
to strongly influence the character of the insulator-metal 
transition. The present study provides a starting point 
for further investigations in doped titanates. 
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FIG. 1. Integrated intensity of the (0.5,0.5,0.5) AF Bragg 
reflection as a function of temperature. The inset shows the 
G-type spin structure. 

FIG. 2. Typical constant-q scans near the AF zone center. 
The upper profile is offset by 200 units. The lines are the 
result of a convolution of a spin wave cross section with the 
instrumental resolution function, as described in the text. 

FIG. 3. Typical constant-energy scans in two different 
directions of reciprocal space. The profiles are labeled by 
the excitation energy in meV. The lines are the result of a 
convolution of a spin wave cross section with the instrumental 
resolution function, as described in the text. 

FIG. 4. Fitted spin wave peak positions in the (1,1,1) di- 
rection of reciprocal space. The line is the magnon dispersion 
curve described in the text. 
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